Wrapping Transition and Wrapping-Mediated Interactions for Discrete Binding along 

an Elastic Filament: An Exact Solution 



David S. Dean,^ Thomas C. Hammant,^ Ronald R. Horgan,^ Ali Naji,^- ^^[^J and Rudolf Podgornik'''^'^ 

^ Universite de Bordeaux and CNRS, Laboratoire Ondes et Matiere d'Aqmtame (LOMA), UMR 5798, F-33400 Talence, France 
''Department of Applied Mathematics and Theoretical Physics, University of Cambridge, 
Centre for Mathematical Sciences, Cambridge CBS OWA, United Kingdom 
''School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran 
^Department of Theoretical Physics, J. Stefan Institute, SI-1000 Ljubljana, Slovenia 
^Department of Physics, Faculty of Mathematics and Physics, 
University of Ljubljana, SI-1000 Ljubljana, Slovenia 
^Department of Physics, University of Massachusetts, Amherst, MA 01003, USA 

The wrapping equilibria of one and two adsorbing cylinders are studied along a semi-flexible 
filament (polymer) due to the interplay between elastic rigidity and short-range adhesive energy 
between the cylinder and the filament. We show that statistical mechanics of the system can 
be solved exactly using a path integral formalism which gives access to the full effect of thermal 
fiuctuations, going thus beyond the usual Gaussian approximations which take into account only 
the contributions from the minimal energy configuration and small fluctuations about this minimal 
energy solution. We obtain the free energy of the wrapping-unwrapping transition of the filament 
around the cylinders as well as the effective interaction between two wrapped cylinders due to 
thermal fluctuations of the elastic filament. A change of entropy due to wrapping of the filament 
around the adsorbing cylinders as they move closer together is identifled as an additional source 
of interactions between them. Such entropic wrapping effects should be distinguished from the 
usual entropic conflguration effects in semi-flexible polymers. Our results may be applicable to the 
problem of adsorption of proteins as well as synthetic nano-particles on semi-flexible polymers such 
as DNA. 

PACS numbers: 87.15.-v (Biomolecules: structure and physical properties), 82.35.Np (Nanoparticles in 
polymers), 82.35.Lr (Physical properties of polymers) 



I. INTRODUCTION 



Interaction, binding and complexation to membranes and polymers, as typified by protein adsorption to DNA 
chains, play a vital role in biology and biochemistry. In biology, for example, a variety of proteins interact with DNA 
in cellular processes involving the reading and packaging of the genome [I] . Transcription factors are proteins which 
bind to specific DNA sequences controlling the flow of genetic information to messenger RNA (mRNA) which in turn 
carries this information to ribosomes for protein synthesis. DNA polymerase is another type of protein (enzyme) that 
plays a key role in the DNA replication process: it synthesizes a new DNA strand by reading genetic information from 
an intact DNA strand that serves as a template. RNA polymerase is a different kind of DNA-binding protein that 
binds to DNA and uses it as a template to synthesize the RNA. In all these examples, protein binding induces local 
deformation of DNA. When multiple proteins bind to DNA they may exhibit a significant degree of cooperativity [2], 
which may result from the formation of loops 3J , specific protein-protein interactions (at short separations) [4] or a 
variable-range cooperative binding of proteins regulated by the tension along the DNA strand [S] . 

Perhaps the most remarkable form of protein-induced deformation of DNA is found in the chromatin fibre. In 
eukaryotic cells, a long strand of DNA is efficiently packed within a micron-size cell nucleus. This occurs through a 
hierarchical folding process [IJ where at the lowest length scale, short segments (about 50 nm) of DNA are wrapped 
around small histone proteins forming nucleosomes: cylindrical wedge-shaped historic octamers of diameter of around 
7 nm and mean height 5.5 nm. DNA-histone binding may be mediated through specific binding sites [5] as well as 
non-specific electrostatic interactions [3i9] since both DNA and histones carry relatively large opposite (net) charges 
(on the average, around 188, 130 and 129 histone charges interact with the wrapped DNA for native chromatin, 
nucleosome core particle and Hl-depleted chromatin, respectively [TU], while DNA itself carries 6 elementary charges 
per nm). In nucleosomes, DNA is strongly bent and wrapped around the core histones in nearly a l-and-3/4 left- 
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handed helical turn, which thus costs a large bending energy given that DNA has an effective persistence length of 
around 50 nm in physiological conditions |14j . 

These DNA-histone complexes are linked together and, on the next level of the hierarchy, fold into a rather dense 
structure known as the 30 nm chromatin fibre, which undergoes a series of higher-order foldings resulting in highly 
condensed chromosomes. Under physiological conditions (salt concentration of around lOOmM NaCl), this fibre 
exhibits a diameter of about 30 nm, while at low salt concentrations the fibre is swollen and displays a beads-on-a- 
string pattern with a diameter of around 10 nm in which core particles are widely separated |^ I151[TB] . In these cases, 
the wrapped structure of DNA around histone proteins remains intact. However, at salt concentrations below 1 mM, 
the DNA begins to partially unwrap from the histone cores [T7j. The electrostatic mechanism for the salt-induced 
wrapping-unwrapping behavior of the DNA around nucleosome core particles has been discussed in the literature for 
general classes of charged polymer- macroion complexes [51 [HI [TS^BU] and for long complex fibres of multiple macroions 
[3T] and will not be considered in this paper. In general, the interplay between the DNA's bending energy cost, the 
electrostatic attraction energy as well as other specific binding energies would be key to understanding the global 
structural features of the chromatin fiber. The precise arrangement of nucleosome core particles in the 30 nm fibre is 
still intensely debated in the literature (see, e.g. [51 IT51 UM f^^HHl] and references therein). 

Genome-wide experimental mappings of nucleosome occupancy in yeast point to a patchy landscape composed 
of at least partly ordered crystal-like configurations with a nucleosome repeat length of about 165 base-pairs [llj . 
It appears that what is encoded in the sequence of DNA is not the sequential ordering of nucleosomes but rather 
the nucleosome excluded regions. In between these regions nucleosomes position themselves via a thermodynamic 
equilibrium mechanism |12| . The ordered configurations of nucleosomes between the excluded regions have been 
modeled as a nonuniform fluid of one-dimensional (ID) hard rods confined by two excluding energy barriers at the 
extremities [TH [13] (see Ref. [ST] for a ground-state treatment in three dimensions) . The interaction between the 
nucleosomes were thus assumed to be of purely hard-core type, i.e., only steric exclusion is taken into account. 
Understanding the details of these interactions and the adsorption-wrapping equilibria for protein-DNA complexes is 
thus of paramount importance. 

Motivated by these findings, we focus here on elastic properties and adsorption-desorption equilibria of a model 
system composed of adsorbing cylindrical particles and an elastic filament modeling a semi-flexible polymer chain 
with short-range adsorption interaction, assumed to be proportional to the arc-length of the filament touching the 
particle. We neglect any electrostatic interactions (which may be justified for charged species at high enough salt 
concentration in the solution) and focus on elasticity and short-range (adhesive) wrapping. Furthermore, we restrict 
ourselves to a 2D Eulerian plane model of a single elastic filament wrapped around one, two or many cylinders. This 
formulation of the model system allows us to introduce an exact formalism based on the Schrodinger representation 
for the partition function of the system and to carry out the calculations explicitly for the wrapping-unwrapping 
behavior of the elastic filament with one or two wrapping cylinders as well as for the effective interaction between two 
cylinders wrapped on a single filament. Our study may thus be applicable to the problem of adsorption of proteins 
as well as synthetic nano-particles to semi-flexible polymers such as DNA. 

In the same context, the tension-mediated interaction between proteins bound to a DNA chain was initially consid- 
ered by Rudnick and Bruinsma [^ within the very same model used in the present paper. The thermal fluctuations 
of the chain were however dealt with only on a harmonic level, corresponding to small Gaussian fluctuations around 
the ground state. Our analysis, based on the Schrodinger representation, goes beyond this approximation and takes 
into account the whole spectrum of chain conformational fluctuations within a rigorous mathematical framework. In 
agreement with Ref. [5], we find that the interaction between two wrapped cylinders may be either attractive or 
repulsive, depending on their relative orientation along the elastic filament. The details of the interaction are however 
different from those found on the simple Gaussian level. 

In another recent related work, Koslover and Spakowitz [35) studied the role of twist in the coupling between bound 
proteins. It was shown that twist resistance results in a more complex interaction between the cylinders exhibiting, 
e.g., damped oscillations superimposed over and counteracting the attractive protein-protein interaction. In this 
paper, confining ourselves to the 2D Eulerian plane, we shall neglect the twist contribution and focus primarily on 
exact solutions for the case where the polymer is described as a twist-free worm-like chain. 

The organization of the paper is as follows: In Section [llj we define our model and discuss the general formalism 
which we shall employ in our study. In Secti on|III[ we apply our formalism to the problem of the unwrapping transition 
for a single wrapped cylinder and in Section |IV[ we study the unwrapping transition of the filament from two cylinders. 
In Section |Vj we consider the problem of effective interaction mediated by the elastic filament's fluctuations between 
two wrapped cylinders. We summarize our results and conclude in Section [Vlj 
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FIG. 1: A schematic drawing of the geometry as well as various quantities defined in the text for one and two cylinders wrapped 
on the Eulerian filament (shown in a symmetric mode). Specifically, for one cylinder (top) the angle at which the filament 
first touches the cylinder (advancing along the arc- length to the cylinder) is ipi and the angle at which it leaves is ip2- For two 
cylinders (bottom), separated by a distance I, the contact angles ai and Q2 are as shown in the figure. 



II. FILAMENT WRAPPING AROUND CYLINDERS: MODEL AND FORMALISM 



A. Filament elastic energy 

Consider an inextensible semi- flexible filament of total arc-length L confined to a 2D Eulerian plane with coordinates 
{x,y). The tangent vector of the filament is normalized to one and can thus be parametrized as t(s) = x(s) = 
(cos(';/'(s)), sin('0(s)). In this parametrization the inextensibility constraint is taken into account exactly. The bending 
energy of the filament is given by 



ds — 
2 



ds2 



ds — 
2 



ds 



(1) 



where k, is the filament stiffness parameter and has dimensions of energy times length (the so-called "persistence 
length" is defined as Lp — K/{kBT) for the ambient temperature T). We assume that one end of the filament is fixed 
at x(s = 0) = and the other end is pulled with an external tension F in the direction x, see Fig. [T] The potential 
energy of the filament is thus given by 



Fx{L) 



-F 



ds cos(7/'(s)). 



(2) 



The total energy of the filament without any interactions with objects in the plane is then given as an elastic energy 
functional Elih] as 



ds 



di}} 
ds 



Fcos(V'(s)) 



(3) 



which is the standard Euler elastic energy expression for the elastic filament under external force. 

The energy function Eq. ([3| also arises in the study of parallel cylinders adhering to elastic membranes |36H39j . 
In this case, when the cylinders are parallel, and if only membrane height fluctuations in the direction normal to 
the cylinders are taken into account, the problem becomes effectively one dimensional and if one uses the Helfrich 
Hamiltonian for the energy of the membrane due to its height fluctuations one can write (up to a constant term) [38j 



E, 



menib 



c 



ds 



dij] 
ds 



a cos(?/'(s)) 



(4) 
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where V'(s) is the tangent angle of the membrane in the direction perpendicular to the cylinders (s being the coordinate 
in this direction), L is the length of membrane in this direction and C the length parallel to the cylinders. The term 
K is the membrane bending rigidity and a the surface tension. The Hamiltonian Eq. Q is thus formally the same as 
that studied by Rudnick and Bruinsma [S] as given by Eq. ([s]); however there are a number of notable and crucial 
differences between the two models. First, Eq. (4]) is multiplied by a macroscopic quantity L, this means that the 
statistical mechanics can be determined by purely energetic (or mean field) considerations by minimizing Eq. Q. 
In this limit one can also assume that fluctuations are small and thus fluctuations can be taken into account via an 
expansion about the minimal energy configuration and treated harmonically as was done by Rudnick and Bruinsma ^S] 
for the polymer model. Another technical difference is that when one considers the interaction between the cylinders 
in the membrane model, the force is calculated as a function of the spatial distance between the cylinders, whereas in 
the polymer problem it is more natural to consider the arc-length along the polymer separating the cylinders as the 
relevant physical parameter. However within the formalism introduced here both ensembles (fixed spatial distance and 
fixed arc-length) can be handled on an equal footing. This choice of ensemble is especially important when analyzing 
the effect of fluctuations. The above points thus conclude the technical motivation for our study: as the problem is 
effectively one dimensional, one can be in a limit where the effects of fluctuations are not merely perturbative and a 
full treatment of fluctuations is thus necessary. The limit where the energetic minimization is valid in the treatment 
of Eq. (3| is when k and F are large (as compared with ksT) while their ratio is kept fixed. 

We also emphasize that our model is one where the size of the objects that are wrapped, the cylinders, is large with 
respect to the microscopic details of the polymer. The effect of absorbed objects on a more discrete model has been 
studied in Ref. [ID], where in particular the effect of absorbed objects on the effective persistence length was taken 
into account. 



B. Wrapping adhesive energy 

Consider now a cylinder embedded in the plane which is free to move and which can become attached to the elastic 
filament. It is important that the cylinder can move since the position at which the filament is attached to the cylinder 
is not constrained in space. We assume that there is a favorable energy of interaction between the cylinder and the 
filament which is local and proportional to the arc-length of the filament touching the cylinder. We further assume 
that in the interaction region between the filament and the cylinder the former follows the surface of the cylinder 
exactly and that the energy of interaction is consequently given by 

El = -jSc, (5) 

where 7 is a line tension parameterizing the affinity of the filament for the cylinder and Sc is the arc-length of the 
filament which is attached to the cylinder. If the angle at which the filament first touches the cylinder (advancing 
along the arc- length to the cylinder) is tpi and the angle at which it leaves is -02 (see Fig. [I] top), then the length of 
the filament touching the cylinder is Sc — R\ip2 — "fpil £^nd thus 

E, = --fR\^2-H- (6) 

Note that if -02 > 4'i the filament wraps anti-clockwise and if ip2 < ipi then it wraps clockwise around the filament. 
As well as having a surface interaction term, there is also a mechanical bending (potential) energy associated with the 
arc-length following the contour of the cylinder. This mechanical energy is given by = E{sc), where the energy 
functional is defined in Eq. ([s]). The path that follows the cylinder arriving at angle 4'i and leaving at angle -02 is 
parametrized over [0, sj as "03 = ("02 — ^i)s/ Sc + ■01- Substituting this trajectory into Eq. ([3]) yields 

= ^1^2 - 01 1 - FR sgn(02 - 01 ) [sin(V'2) - sin(V;i)] . (7) 

The total energy of the system is then composed of a bulk elastic term, Eq. ([S]) , for the free segments of the filament 
and a boundary or contact interaction energy Ec with the cylinder which is a sum of the terms in Eqs. (|6]) and ([7|, 
i.e., Ec = Ei + Em- 
One notes that the above energy functional is analogous to a ID Coulomb fluid with a charged boundary, but with 
an imaginary surface charge [HI |32 . 



C. Partition function of a filament-cylinder complex 



To determine the conformational equilibrium of a system consisting of an elastic filament with any number of 
wrapped cylinders we must evaluate the corresponding partition function. To do this we first define the evolution 
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kernel K {tfj , tfj' , s) which evolves the state of the filament for a distance s measured along its arc-length with the 
boundary condition that '0(0) = '0 and ip{s) = V''- This kernel is defined by 

K{^,ij',s)= / d[i:]exp{-l3E[i,]) (8) 

JV(o)=i/' 

where (3 = and £'[■0] is the elastic energy functional in Eq. (|3|. It determines the evolution in arc-length s of 

a wave Junction /(0) by 



d0'if (V', 0', s)fm = exp(-sF)/(V'), (9) 
where H is the corresponding Hamiltonian operator 

The partition function for a single cylinder, which first becomes attached to the filament at arc-length s — li measured 
from the fixed end x(s = 0) = (see Fig. [T]), is then given by 



Z = j (i'0od'0lrf02rf'03-?^('0O,'01,^l)5'[01,02]-?f(02,03,-^ - - -R|02 - -011), (H) 

where 0'(O) = ipQ and '4'{L) — ip^ and the definition of 5* follows from the contact energy (or the interaction energy 
between the cylinder and the filament as defined in the preceding Section) and is given by 

5[^i,V2] =exp(/3|V'2-0i|(7^- ^) + /3i^^sgn(0;2-V'i)[sin(02)-sin(V'i)]) . (12) 

The generalization to a system in which the elastic filament interacts with several wrapped cylinders is obvious from 
the above. 

In the limit of a long elastic filament where the cylinders are close to its midpoint, we make use of the fact that the 
propagator at large arc-lengths is given in terms of the ground-state wave function of the corresponding Schrodinger 
equation 

i^(V, V'', L) = exp(-Soi)^'o(V')*o(0''), (13) 

where Eq is the ground-state energy of the Hamiltonian H, Eq. ( |10[ ), and the corresponding wave function, i.e. 
H'^Q^ip) = Eo'^q{iI}). The ground-state energy in the ground-state-dominance limit then gives also the partition 
function in the limit of large arc-lengths. 



D. Constrained and free wrapping of cylinders 



Because of specific interactions between the filament and the cylinders the wrapping angles "01 and 02 may not 
be completely free and constraints may exist that reduce their overall degrees of freedom. For instance, in the limit 
where the adhesion energy 7 is very large, the elastic filament could wrap around the cylinder a number of times, via 
an escape into the third dimension that does not cost any extra energy, or so we assume in this study. 

The simplest constraint is to take the wrapping angle (see Fig. [l] top) "02 — V'l — ol\ fixed. The constrained 
partition function at a fixed wrapping angle can be computed and then the ensemble where it varies can be obtained 
by carrying out the remaining integrations with the appropriate statistical weights. This holds for a single wrapped 
cylinder but can be extended to the case of several cylinders as well. 

For a single cylinder we then define the constrained partition function for the ensemble where "02 ~ V'l — cti is kept 
fixed as 



Z{ai) = J dtljodipidtlj3K{tljo, 0i, Zi) exp {/3FR sgn(ai) [sin(V'i + ai) - sin(0i)] ) x 

X K{ipi + ai,i)3,L- h - R\ai\). 



(14) 
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The partition function for the ensemble where the contact angle can vary freely is consequently given by 

rL-h/R 
-L-h/R 



z = 



daiZ(ai)exp(/3|ai|(7i?- ^)) . (15) 



Note that given that the point of first contact between the cylinder and filament is at s = ^i, the maximum amount of 
filament that can be wrapped about the cylinder is L — Zi, corresponding to wrapping angles of ±(i — (clockwise 
and anti-clockwise). When ai is fixed, the surface terms are periodic in ipi and one can use standard Mathieu function 
analysis to carry out the computations. The procedure is described in detail in Appendij|A] 

For two cylinders separated by a distance Z, where the first one is located an arc-length li away from the fixed end 
and with contact angles ai and a2 (see Fig. [T] bottom), we have in complete analogy with above 



Z(ai,a2) = J dt/^odipid^padip^ K{ipo,^i,li)SR^{ipi,ai)K{'il;i + ai,ip3,l) X 

X {ip3,a2)K{ip3 + a2, ips, -^12), (16) 

where 

L12 = L - I - h - Bi\ai\ - R2\a2\ (17) 
is the arc-length between the right, free, end of the chain and the rightmost cylinder (see Fig. [l] bottom), and 

SR{'tp, a) = exp {jSFR sgn(Q!) [sin(V' -I- a) — sin(V')]) , (18) 
where i?i is the radius of cylinder 1 and R2 is that of cylinder 2. 

E. Constrained arc-length separation between wrapped cylinders 

Now consider the limit of a long filament where the cylinders are close to the midpoint, i.e., we take li and L — li 
large while keeping the arc-length between the final point of the filament touching cylinder 1 and the first touching 
cylinder 2, Z, constant (see Fig. [l] bottom). These are all lengths along the chain (arc-length) which are fixed, if we 
wish to keep the physical distance in a given direction fixed we must consider another ensemble; we shall discuss this 



later. In this limit (corresponding to the limit of ground state dominance, Eq. (13)), we find, up to an overall factor. 



Z[ai,a2) = exp(-So-Z^i2) j dipid^3^Q[tpi)SR^{'4)i,ai)K{tpi -I- ai, -03, O'^i?,^ ("03, a2)*o(V'3 +a2)- (19) 

An interesting case emerges when the two wrapping angles are equal and correspond to complete single wrapping, 
OLi — a2 = 2tt. The I dependent part of the partition function is then given by 



Z(0 = cxpiEol) I dVirfV'3*o(V'i)^(0i + 2^, ^3, 0*0(^-3 + 2tt), (20) 
and using the periodicity this gives 

Z{1) = exp(^oO / d0id03«'o(^i)if(^i, ^3, 0^0(^-3). (21) 



The same result is found if one assumes that the wrapping is antisymmetrical, a2 — — ai = — 27r. The symmetry 
of wrapping between the two cylinders will be addressed furthermore below. The difference between symmetric and 
antisymmetric wrapping is obvious from Figs. [l]and[2] 



F. Wrapping-unwrapping transition of cylinders on an elastic filament 

We first consider the case of a single adsorbing cylinder. The elastic filament is of length L and a wrapped cylinder 
is positioned at li. We restrict ourselves to the case where both L and li — > 00, assuming that li is not fixed so 
that the cylinder can attach anywhere and consider the case where the elastic filament can wrap on the cylinder any 
number of times. 




FIG. 2: Top: Schematic view of the wrapping of two cylinders in an antisymmetric mode where the wrapping angles for the 
two cyhnders differ in sign. Bottom: the looped configuration (full line) of two wrapped cylinders corresponds to a negative 
value of the horizontal projected separation d± and the extended configuration (dotted line) with positive d±. 



In this case, the expression in Eq. (14) becomes, up to a constant prefactor, independent of a: 



Z{a)^exp{^Eo{L^R\a\))f{a), with f{a)= / ^'2(V')5fl(V, «), 

Jo 



(22) 



and the function f{a) is clearly bounded. From Eq. ( |15[ ) the partition function for the ensemble with variable a is 
obtained by integrating Eq. (22) over a from —L/R to +L/R^ i.e. 



/' / K RE 

da exp /3|a|(7i? -7^ + — ^) ) f{a). 
-L/R \ 



2R 



Defining 



AE = jR- 



K REo 
2R ^ 



(23) 



(24) 



we see that the elastic filament will then wrap around the cylinder a macroscopic number of times (in the sense 
that the total length wrapped around the filament will be of the order of the filament length) if AE > 0. However, 
for AE < the filament will have only a microscopic length wrapped around the cylinder. The equation AE = 
therefore defines a wrapping transition in the phase diagram of variables (3 and as a consequence of competition 
between the wrapping energy of the cylinder and the configurational entropy of the chain. 

In the zero-temperature limit where /3 — ^ 00 the ground state of the filament configuration is a straight line in the 
direction of the applied force and so without a cylinder 



Z{L) w cxp{-LEq) = exp{l3FL), 



(25) 



or Eo{T = 0) = -~/3F. Therefore, if -fR — k/2R — FR > 0, the system is in the wrapped state, otherwise the filament 
unwraps from the cylinder. This conclusion is easy to deduce from purely energetic arguments for large values of the 
wrapping angle a. 

We next consider several cylinders in the plane and assume that the filament can wrap around all of them without 
impediment. The wrapping will induce effective interactions between the cylinders that can be either repulsive, 
attractive or non-monotonic. Wrapping transitions with effective interactions between cylinders can be viewed as a 
model for nucleosomal wrapping [T21[I3ll^- However, the details of the most general case remain to be elaborated. 

For clarity we briefly describe the system with three cylinders: cylinder 1 at position ^qi from the left end of the 
chain, cylinder 2 separated from 1 by a distance I12, cylinder 3 separated from 2 by a distance I23 and with contact 
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angles ai, ai and 03. The partition function can be written as 

Z(ai,a.2,a-i) = J d^odil}idip3d^5dip7K{ipo,iJi,loi)S{iJi, ai)K{'ipi + ai, i/'s, ^12)5(V'3, "2) x 
xK{tp3 + a2,tp5,h3)S{ip5,a3)K{tp5 + a3,tpr,L- {ki + I12 + hs, + Ri\ai \ + i?2|a2| + ^slasD), 

(26) 

with the general definition 

'5(V'2i+i, tti) = exp {pPRi sgn(ai) [sin(V'2i+i + a^) - sin('(/'2i+i)]) , (27) 

where i (for i = 0, 1, 2) is the index of the wrapped cylinder. 

The general system for N particles bears some overall similarity to the ID Tonks gas, although the effective 
interactions between the cylinders are more complicated. In order to make this system applicable to the problem of 
nucleosome wrapping around DNA, it may be necessary to insert a chemical potential term for the cylinders similar 
to the case of a grand-canonical Tonks gas [T^lIIS]- The main difference between the two models is, however, in the 
fact that the interaction between the wrapped cylinders along the elastic filament can be much more complicated than 
in the Tonks case and could in principle lead to "phase separation" without imposing any excluding energy barriers 
along the chain |21) . 



III. UNWRAPPING TRANSITION: ONE CYLINDER 



We now apply the general theory derived in the previous Sections to the problem of the unwrapping transition for 
a single wrapped cylinder. 

To accommodate the wide range of values that the variables take in the physical systems to which our theory would 
be applicable, it is convenient to recast the calculation in terms of equivalent dimensionless variables. We consider an 
elastic filament of persistence length Lp, fixed at one end and pulled from the other end with a force F. The wrapped 
cylinders are all of radius R. We define the following dimensionless variables 

or 

11^ —f, f ^ PRF, cj^ P-IR, e = ER, (28) 

where Lp — /3k and 7, the wrapping energy of a cylinder per unit of arc-length, and other parameters have been 
defined before, Eqs. ([T])-([5]). In terms of these variables the Hamiltonian can be written a.s H = H' /{jiR) and thus 
H' has eigenvalues which can be written in the form of the Schrdinger equation 

i/Vm = MemV'm (29) 

where = E^R. 

As an illuminating example let us consider a specific case with the dimensionless quantities /x = 1 and a = 0.5. 
Assuming cylinder radius of i? = 2 nm, this gives actual parameter values as Lp = 1 nm (persistence length), 
K = ksT X ~ 4 X 10~'^°J.m (bending stiffness) at T = 300 K, and 7=1 pN (line tension for cylinder-filament 
adhesion) . 

For the single cylinder wrapping transition, a rough estimate for its occurrence can be obtained in the following 



way. In Section II F it was argued that for large /? (low temperature) Eq ^ fiF and that the wrapped phases occurs 



for forces which obey /S.E > 0, Eq. [24j For scaled quantities as introduced above this inequality becomes 

/ < cr - ^ and thus ct' = cr - ^ > . (30) 



Alternatively, we can then consider the Schrdinger equation, Eq. 29 (see also Eq. (|10|), as an approximate oscillator 
equation and get an improved estimate Eq = — /3F + \J F j (4k). Using Eq. po| above, the condition for the 
wrapping transition then becomes 



/-y772^-^'<0, (31) 
This can be recast as a quadratic equation in -^/J and the condition can be satisfied if the discriminant is positive, i.e. 

ct' > - — and thus a > - - — . (32) 
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FIG. 3: The wrapping transition of a single cylinder on an elastic filament. Exact solution for the average wrapping angle ratio 
<OL> jamax is shown as a function of the dimensionless external tension / = pRF for fixed /i = 1 and a = 0.75. 



In the limit where is large, the above expression reduces to the condition in Eq. (30). 
corrections to this form due to the non- harmonic terms in the effective potential in Eq. ( 29 ) 



There will of course be 
Nevertheless, from both 

arguments it follows that a sufficient condition for the transition to occur is ct > /^/4. If we return to original variables 
we see that the second condition of Eq. ( 32 ) gives 



7 > 7 



where 



7 = 



1 K 

2m 



16k 



(33) 



In the limit of zero temperature this result agrees with that of Weikl [36J for a cylinder and membrane. The interpre- 
tation of the zero-temperature result is simply that the adhesive energy must overcome the bending energy to make 
the bound state stable. We see that the effect of temperature is to diminish the minimum adhesion energy necessary 
to obtain a bound phase. At first sight this appears counter-intuitive and so we investigate the effect of temperature 
below in more detail. 

We first present an exact calculation of the wrapping transition in a particular case. The mean value of a, i.e. (a), 
is calculated for fixed {n,a) as a function of the reduced force /. Whilst taking the filament to be of infinite extent 
we restrict the range of a to be —a^ax < a < OLmax and take amax = 100. For small / we expect (a) to be large, 
which corresponds to the elastic filament maximally wrapped onto the cylinder: \a\ = amax- As / is increased the 
system will go through a transition from the wrapped to the unwrapped configuration. 

To solve the corresponding Schrdinger equation, the range of x is chosen to be an integer multiple N of 27r, where 
N would then be the number of lattice sites in the Bravais lattice for the discretized Hamiltonian. The Schrodinger 
equation (Eq. |29[ ) is then recast as a matrix equation with imposed periodic boundary conditions. This will lead 
to a band structure with the Brillouin zone and the number of Bloch states per band determined by N. Since the 
solution of this problem requires only strictly periodic eigenfunctions we choose A'^ = 1 and the range of x as 2tt. The 
eigenvalues are then equal to iiem- 

In Fig. [3] we show the results for the mean wrapping angle <a>, scaled by amax, as a function of the external 
dimensionless force / for fixed /i = 1 and a = 0.75. The wrapping transition, corresponding to the maximum in the 
derivative d{a)/df, then occurs for the critical reduced force fc — 1.2 which, for the parameters stated above (i.e., 
using Lp = 1 nm, R = 2 nm) gives Fc ~ 2.4 pN. Compare this with the corresponding estimates from 



(i) Eq. (30). We find 



fc 



Ai/4 



0.5 



Fc = 1 pN. 



(ii) Eq. (31 ). We solve the quadratic equation (for /i = 1 and a = 0.75) to get 

fc = 1.31 Fe = 2.67pN, 
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FIG. 4: The wrapping transition as a function of the actual force F and for temperatures T = 210 K up to 390 K in steps 
of 30 K. The physical values for the parameters are given in the text. In both cases the curves move right as T increases. 
Left pane: we have a / fj, = 0.15, a' < 0. The wrapping transition is relatively smooth and as T increases it becomes stronger 
and moves to larger values of F. Right pane: cr//i = 0.35, a' > 0. The wrapping transition is much stronger and sharper as 
compared with the former case and, as predicted, the transition moves to larger values of F as T increases. It is this behavior 
in both cases which might be considered counter-intuitive. 



Obviously, in order to get an accurate prediction for the critical force it is necessary to solve the full Schrdinger 
equation as we did above. Nevertheless, the harmonic approximation in Eq. (31 ) is in acceptable agreement with this 
exact result. 

In order to highlight the effect of temperature it is important to consider the dependence of (a) on the force F 
rather than on the dimensionless reduced force / because the latter contains a hidden dependence on T which would 
obscure the effect we are studying. In the harmonic approximation, we have that the critical force F^. is predicted to 
be 



'Fr 



4^^ 



+ (7 - 7') 



(34) 



We investigate two cases and take typical values for the parameters similar to those stated at the beginning of this 
Section, i.e., R = 2 nm, Lp — 1 nm, k = ksT x Lp and the temperature is taken as T = 300 K, so that n — 1. We 
study two cases corresponding to a/fi — 0.15 and 0.35, respectively. We also note that these choices respectively give 
a' > and a' < 0, which allows us to test the significance of the inequality in Eq. (32 1. We study a wide range of 
temperatures above and below T = 300 K which is necessary to reveal the T-dependence of the wrapping transition. 
Of course, such a wide temperature range does not occur in vivo but our aim is to test the prediction inferred from 
Eq. ( [M] ), namely that the critical force for the wrapping transition increases, rather than decreases, with increasing 
T. The values for T considered are between T — 210 and 390 K in intervals of 30 K, and the results are shown in 
Fig. |4]for the above two cases: 

(i) a/fj. — 0.15, a' < 0: This corresponds to relatively weak adhesion coefficient 7 = 0.30 pN and whilst we expect a 
wrapping transition to occur, we expect it to be relatively smooth and, according to Eqs. (33 1 and (34), to turn 
on more strongly as T increases and for the critical force to increase as T increases. We see these features in 
Fig. |4j left pane. We conclude also that even though a' < the wrapping transition does occur. For T — 300 K 
the harmonic approximation in Eq. (34) predicts Fc — 0.52 pN. The transition is not sharp and so it is not 



feasible to read a specific value for Fc from the graph but it is clear that this prediction is very much on the low 
side. 

(ii) a/n ~ 0.35, a' > 0: The adhesion coefficient 7 — 0.70 pN is larger and, in this case, we expect the transition 
to be stronger and sharper than in case (i). This is seen in Fig. |4j right pane. As T increases the transition 
moves to larger critical force as predicted. For T = 300 K the harmonic approximation in Eq. (34) predicts 
Fc = 1.37 pN. From the graph the value of Fc follows as f c ~ 3 pN, twice as big as the harmonic prediction. 
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We have thus verified that the effect of increasing temperature T is to increase the critical force at which the wrapping 
transition occurs and this feature is clearly seen in Fig. |4]jii)- In other words we can conclude that increasing the 
temperature causes the wrapping transition to occur at larger force and the corollary is that for fixed force increasing 
the temperature can actually cause the transition: a somewhat surprising result! 

These effects can be understood by noting that it is the ground-state energy for the full Schrdinger equation 



which contributes in Eq. (22 1 and this is clearly T-dependent and different from i?o for the harmonic approximation. 



As T increases the potential becomes narrower leading to Eq increasing with T. This leads to a smaller threshold, 7', 



as defined in Eq. ( 33 ) 



IV. UNWRAPPING TRANSITION: TWO CYLINDERS 

We now study the unwrapping transition in a system of two cylinders of radius R, labelled by i = 1, 2, wrapped by 
a length of filament. There are two cases which we consider. 

The first case (case I) consists of the two cylinders wrapped by a fixed length of filament. The relevant variables, 
shown in Fig. [l] are as follows: 

- The total angle, a^, wrapped by filament for the i-th cylinder. These are both fixed (quenched) for i — 1,2. 

- The length, I, of filament directly between the two cylinders, i.e., from the exit of the first cylinder to the 
entrance of the second, which is also fixed. 

We encode the configuration of the cylinders by measuring the expectation value of d_\_, the projected horizontal 
displacement between the centers of the two cylinders. The maximum projected horizontal displacement is given by 
dmax = I + 2i?. We shall deal with this case in Section |IV A| below. 

The second case (case II) consists of the two cylinders pinned a distance /' apart along a filament, which are wrapped 
by a dynamical (or annealed) length of filament. The relevant variables are: 

- The length V of filament between the sites of pinning of the two cylinders. This is fixed (quenched). 

- The total angle, a^, wrapped by filament for the i-th cylinder. These are dynamical (annealed) and free to 
change. 

- The angle, a^, wrapped by the internal length of filament (initially of length V) between cylinders for the z-th 
cylinder. These are also dynamical; they take into account the ways in which the filament can wind. 

- The length, Z, of filament directly between the two cylinders, i.e., from the exit of the first cylinder to the 
entrance of the second. It is now dynamical being determined hy I = I' — R{a[ + a2)- 

Again we encode the configuration of the cylinders by measuring the expectation value of dj_, the projected horizontal 
displacement between the centers of the two cylinders. The maximum projected horizontal displacement is given now 



by dmax ^ I' + 2i?. We shall deal with this case in Section IV B below 



A. Two cylinders: constrained wrapping 

In the first case (case I) where the wrapping angles are fixed, we write the constrained partition function of this 
system as 

Z{ai,a2,l) = y"dVirfV'2(0|O(V'i,V'i+ai)exp(-iI//i?Ai)O(V'2,V'2 + a2)|0), (35) 

where we assume that the ground state dominates the external regions of the elastic filament, i.e., regions outside 
the part bounded by the two wrapped cylinders. The operator insertion O corresponds to the wrapping of the elastic 
filament around a cylinder and is given by 

0(V,,V + a) = C„(V)|V)(V' + a|, (36) 

where 

CaW = exp (|a| [cr + eo] + / sgn(a) [sin {tp + a) - sinip]) . (37) 



This gives for the constrained partition function Eq. ( 35 1 
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FIG. 5: Dependence of the average horizontal separation (dx) (normahzed to dmax ~ I + 2R) on the external dimensionless 
tension / — jBRF for different sets of parameters in the case of constrained symmetric wrapping with ai — a2 — it. Left pane: 
fi = 10 (high rigidity). Right pane: /.i = 5 (low rigidity). Negative values of the average horizontal separation indicate the 
presence of a looped phase, i.e., the cylinder at larger distance along the elastic filament lies to the left of the other cylinder. 
Curves from left to right correspond to / = nR for integer n — 3, ... ,8. 



Z{ai,a2,l) = / dV'idV2^'o(V'2)C„i(^i)^e-^"'/^P™(V'i+«i,V2)C„,(V^2)*o(^2 + a2), (38) 

where 

Pm{^l,^2) = «'m(V'l)*r„(V'2). (39) 

The functions Pm and the exponential in front of the Pm can be pre-computed outside any nested integration loops for 
fixed ai, leading to a quick and straightforward evaluation of the above expression. We are interested in the average 
horizontal distance between the centers of the two cylinders which is given by 

{d±) = — ^ /ds /d(sin^,)/d^id^2f(0|O(^i,i^i-f ai)e-^^/^'^|i^,) X 

Z(ai,a2,l) J J J \ 

X (V'.|e^^('-^)/^^d(V'2>2 + a2)|0) +i?sgn(ai)(sin(i^2)) -i?sgn(a2)(sm(V'3)) 

(40) 

A form suitable for numerical computations is reworked in Appendix [B| 

We remark at this point that in our model nothing stops the cylinders from passing through each other. Therefore 
the inter-cylinder force for small horizontal separations has to be interpreted only up to the limit of cylinders actually 
touching. The definition of this point depends on the symmetry of the wrapping, i.e., it differs when ai — a-i and 
when a\ — —ai. These arc the two cases that we have designated symmetric and antisymmetric, respectively. 

The wrapping mediated interactions derived in this way should then be added to the hard core ID Tonks gas model 
[43j , which is what was only taken into account in the studies of the positional distribution of nucleosomes along the 
genome [TH[T3]. We do not delve into this problem specifically here, but plan to address it elsewhere. 

An interesting case with two cylinders was discussed in detail by Rudnick and Bruinsma [S]: they dealt with a 
system composed of two cylinders at a fixed arc-length separation / and with fixed wrapping angles, solved in the 
Gaussian approximation level. We analyze here the exact solution for the two phases described in Ref. ^ which are 
(see also Fig. [2]): 

(i) the looped phase where the mean projected separation (dj^) is negative meaning that the cylinder at larger 
distance along the filament lies to the left of the other cylinder, thus causing the elastic filament to loop and 
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FIG. 6: Dependence of the average horizontal separation {d±) (normalized to dmax ~ I + 2R) on the external dimensionless 
force / = PRF for different sets of parameters in the case of constrained symmetric wrapping ai — 02 — Svr/S. Left pane: 
= 10 (high rigidity). Right pane: jj, — 5 (low rigidity). Clearly, for this set of parameters only the extended phase is allowed 
for large values of force regardless of rigidity. Curves from left to right correspond to I — nR for integer n = 3, . . . , 8. 



(ii) the extended phase where {dj_) is positive and there is no loop. 

We fix Qfi and a2, and choose a range of elastic filament arc-lengths, I, between the cylinders. We calculate the average 
horizontal separation {d±) for a given external tension / (in dimensionless units) and plot / versus {d±) /{I + 2R) for 



various choices of parameters. The average horizontal separation is calculated from Eq. (40). 

For symmetric wrapping described by ai = a2 = tt and /i = 10 (high rigidity) the results are shown in Fig. [s] left 
pane. From left to right the lines are for I — nR for integer n — 3 . . .8. For small forces the filament is in the looped 
phase. Increasing the force actually makes the loops larger as shown by {d±) becoming more negative. Eventually 
the extended phase becomes preferable and we have a limit of full extension as the force increases. Small I also makes 
the loops more energetically favorable. For ai = a2 = tJ" but smaller (rescaled) rigidity, /x = 5 (Fig. [s] right pane), 
we see that only the extended phase is allowed for the largest values of I. This implies that for a given rigidity there 
is a minimum distance (and obviously maximum force) for loop formation. 

The largest extension in both these cases is noticeably less than the hypothetical maximum dmax — I + 2i?. This 
can be understood for the chosen fixed wrapping angle of tt since high rigidity constrains the filament to be close to 
tangential at both the entry and exit points. 

In Fig. |6]we then plot / versus {d±)/dmax where dmax = I + 2i? for 11 — 10 (left) and fi — 5 (right) with fixed 
Q!i = Q!2 = 37r/8 < 7r/2. The curves for the two rigidities are qualitatively similar with only the extended phase 
allowed. There is a small / dependence, with curves corresponding to shorter I lying to the right of those for longer 
I. Unlike in Fig. [5] the largest extension for both values of /i is now closer to dmax since the fixed wrapping angle is 
smaller than tt and the constraint that the filament is tangential to the cylinders has a weaker effect. We observe that 
as the rigidity decreases the dimensionless displacement for a given (small) external tension also decreases; a lower 
rigidity allows the elastic filament to fluctuate more and increase the entropic contribution to the free energy without 
incurring a large energy penalty. From our exact analysis, we thus qualitatively confirm the results for this system by 
Rudnick and Bruinsma |S]. 

We will not discuss in detail the antisymmetric configuration 02 — —ai as here a looped phase does not occur and 
the plots are similar to those shown in Fig. [6] 



B. Two cylinders: unconstrained wrapping 



We now allow the to be a dynamical or annealed variable as opposed to the quenched case of the previous 
Section. As noted in the beginning of Section IV (case II) , the corresponding realization would be an elastic filament 
wrapped around two cylinders with a fixed arc-length separation I' between the two points of pinning and no other 
constraints. The length of elastic filament adhering to the cylinders is now variable and so its length outside the 
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FIG. 7: Mean projected separation, {d±} (blue curves), and mean wrapping angle, (a) (red curves), for unconstrained wrapping 
of elastic filament around two cylinders. Left pane: a = 1. Right pane: a = 10. Solid curves correspond to = 10 and dashed 
curves correspond to /x = 1. Note that the results have been normalized to maximum displacement, d^ax = I' + 2R, or to 
maximum angle. 



region between the two cylinders and the effective length between them must both be allowed to vary dynamically. 

In Section IV A we gave an expression for the partition function, Z{ai, a2,l), for fixed angles of wrapping and 
fixed effective length of elastic filament between the cylinders. Therefore the total partition function for dynamical 
wrapping is clearly given by 



min(Z/i?,Qi ) pmin(l/ R—a'^.a2) 

daida2 I da'^ I dai^ Z{ai,a2,l = I' — Ra'^ — Ra2), (41) 



where the integration is over the four wrapping angles which take into account the various ways for the wrapping to 
occur; for the i-th cylinder ai (or a^) is the total (or internal) wrapping angle of the elastic filament, where internal 
refers to the length between the two cylinders. The computations were done with ai,a^ > 0; although the choice of 
sign for the first cylinder wrapping is arbitrary, we force symmetric wrapping on the second cylinder. We find that 
the antisymmetric configuration gives similar results and so do not report on it in detail. {d±) can now be calculated 



by averaging the expression in Eq. (40) over the four wrapping angles. 

Because of the ability to pre-computc the various contributions to the partition function before any integration 
and/or summations are done, computation time is kept to a minimum. However, there is an increase of several orders 
of magnitude in the computing cost compared with the fixed wrapping case. 

In the following computations, we chose I' = 27ri? and < 127r. This puts a reasonable limit on the maximum 
wrapping angles but it already has a significant outcome and the amount of computing resources required is kept 
reasonable. 

In Fig. [t] we plot the average projected length {d±) versus / for cr = 1, 10 and n = 1, 10. As / is increased an 
unwrapping transition is indicated by a rapid decrease in (a) and a corresponding increase in {dj_). The results for 
a = 1 show a clear unwrapping transition (Fig. [7j left). When the unwrapping occurs, the distance between the 
cylinders increases as one would expect, whilst in the wrapped phase there appears to be a small but potentially 
interesting effect causing a decrease in the separation as the force increases. This can be seen by the bump in the 
fi = 10 line at small /. As the rigidity fi is decreased the unwrapping is less pronounced; it is now easier for the 
filament between the cylinders to fluctuate more strongly on average. 

For the case of cr = 10 (Fig. [7j right), there is no unwrapping transition for both values of fj, over the range of 
forces indicated. One would expect this kind of outcome when one changes the binding energy a. Higher binding 
energy prevents unwrapping, holding the cylinders tightly onto the filament. The distance between the cylinders does 
however slowly increase with increasing rigidity owing to the straightening of the elastic filament. 

We also note, in this case, that, although a does not change (it stays maximally wound), the average projected 
length between the two cylinders decreases for large enough external force. This counters the fact the filament will 
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FIG. 8: Free energy as a function of the dimensionless external tension / for fi = 10 and a — 4.5. Left: Red solid curve 
corresponds to a single cylinder a — tt, dashed curves correspond to antisymmetric double cylinder, oi = — a2 = tt. The 
lower (blue) curve corresponds to / = 2nR and the upper (black) curve is limit of increased separation. Right: Red solid curve 
corresponds to a single cylinder a — tt, dashed curves correspond to symmetric double cylinder, ai = 02 = n. The upper (blue) 
curve corresponds to I = 2nR and the lower (black) curve is limit of increased separation. Note the double cylinder (black) 
curve that coincides with the single cylinder (red) curve at free energy equal to zero gives the limit of infinite separation. 



tend to straighten with increased external tension. It does, however, lead to the conclusion that for large external 
tension there is an effective attraction between the cylinders that pulls them together. This conclusion is corroborated 
also by direct evaluation of the interaction free energy between the cylinders as discussed below. 

C. Free energy of the unwrapping transition 

To understand the phase transition between a free and wrapped cylinder, i.e., the phenomenon of unwrapping and 
desorption, we calculate the free energy of wrapped and unwrapped systems. We restrict ourselves to the case of fixed 
wrapping angle (case I). This can be trivially calculated from our theory as we already have calculated the partition 
fimction: 

P<t>{ai,a2,l) = -\nZ{ai,a2,l). (42) 

Since we normalize all energies by subtracting the ground-state energy and do not put in wavcfunctions for the end 
points, we normalize the system to have Z = 1, ^ = for a filament with no wrapped cylinders. We consider the free 
energy of a system for the cases of a single wrapped cylinder with wrapping angle fixed to a = tt, and of the system 
of two wrapped cylinders in the symmetric and antisymmetric configurations with ai = n and a2 — ±7r, respectively. 
In Fig. |8] we show the free energy for these three cases. The lefthand and righthand plots compare the single cylinder 
with the double cylinder for the inter-cylinder separation / = 2'kR and the limit of large I for the symmetric and 
antisymmetric configurations, respectively. We use /i = 10.0 and a = 4.5. For / = 0, the free energy of the two 
cylinder system is essentially double that of a single cylinder, as would be expected. As the force increases so does 
the free energy. 

In the antisymmetric case, at low external force the wrapping of both cylinders is energetically favorable. As the 
external tension increases the two cylinders stay wrapped until both desorb simultaneously (the line crosses the x- 
axis). As we would expect in the limit of large separation, the critical force for unwrapping is the same as that of a 
single cylinder, however as the separation decreases a larger force is required as they are bound together. 

In the symmetric case we also see that at low external tension the double cylinder wrapping is energetically favorable. 
However, as the tension increases there comes a point where the single cylinder wrapping becomes energetically 
preferred and one of the cylinders then unwraps and leaves the chain. The details of this process can not be captured 
appropriately in our model since we do not include the chemical potential for the unwrapped cylinders. 

Increasing the external tension further unwraps the last cylinder. If the two cylinders are far apart we see the same 
behavior as for the antisymmetric case, as would indeed be expected. This agrees with previous studies showing that 
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antisymmetric (symmetric) wrapping gives an attractive (repulsive) force [5]. 

V. EFFECTIVE INTERACTION BETWEEN TWO CYLINDERS 

We now consider the problem of effective interaction mediated by the elastic filament's fluctuations between two 
wrapped cylinders. The elastic and adhesive energy expressions are given as before. 

A. Constrained wrapping angles 

We first consider the case where the two cylinders are separated by fixed arc- length, and the wrapping angles 



are constrained at fixed values ai and a2 (see the remarks identifying case 1 in Section IV and IV A I. In numerical 
simulations we will assume furthermore that they are equal up to the sign, ai = ±q;2; the sign differentiates between 
the symmetric and the antisymmetric wrapping case. We can evaluate an effective interaction between the cylinders 
by fixing the (projected) horizontal separation d_L between them and then calculating the corresponding free energy. 
The projected separation in the embedding space is defined as 

d_L — i? sgn(Q;i) sin('0i + ai) + / ds cos('0(s)) — i? sgn(a2) sin(?/;2)- (43) 

This constraint can be handled most conveniently by introducing an additional term into the total energy of the system 
via a Lagrange multiplier of the form Xd± , where A can be interpreted as the force used to impose the constraint of 
fixed projected distance between the two cylinders in the horizontal direction. 

The wrapping constraints can be implemented as before and one finds the appropriate partition function to be 

Zx{ai,a2,l) = y"dVirfV'2(0|OA(^i,V'i+ai)e-^^'/^^OA(V'2,V'2+a2)|0), (44) 

where 

Hx'i'Xrniij) = (-4t^ - K.f - A) cos (V')) *Am(V') = M^Am^'Am (V') , (45) 



with 



and 



Oa('0i,V'i +ai) = O('0i,'0i +Q:i)exp(-/3Ai?sgn(ai)sin(-0i +ai)), (46) 



Oa(V'2, iI^2 + a2) = 0(V'2, V'2 + a2) exp (^Ai?sgn(a;2) sin(-(/'2)) ■ (47) 
With this constraint the partition function Z{di_) becomes 

Z(d^,«i,a2,0= / — e^^^-ZA(ai,a2,0- (48) 

The appropriate constrained free energy then follows as 

/3ri(rf_L,ai,a2,0 = -lnZ((ij_,ai,a2,0- (49) 

The generalization to many wrapped cylinders is formally straightforward but computationally very tedious. One 
interesting future endeavour would be to assess the effects of non-pairwise additivity in the case of constrained and 
unconstrained wrapping around the interacting cylinders, i.e., the dependence of effective two-cylinder interaction on 
the presence of other wrapped cylinders along the elastic filament. 

We can define related free energies appropriate for the sake of numerical calculation using 



/?$(«!, a2,^, A) = -lnZA(Q:i,Q;2,0 



(50) 
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FIG. 9: Free energy as a function of the horizontal displacement projection (dx) is shown for the case of constrained wrapping 
with oi = ±0:2 = 57r/8, / = 0.4, 1, 2, 3 (from lower to upper curve), = 50 and a = 13.0. In the antisymmetric case (left) the 
effective interaction is attractive and in the symmetric case (right) it is repulsive. 



and its Legendre transformation 

S(;, (d_L),Q!i,a2) = min($(ai,a2,^, A) - A(d_L)) , (51) 



A 



which correspond to a system where either the polymer length, Z, or the average horizontal displacement, are 
fixed respectively. 

In Fig. |9]we plot $ for A = 0, ai = Stt/S and 012 = ±Q;i corresponding to symmetric and antisymmetric config- 
tirations, respectively. The length of elastic filament connecting the cylinders, I, is not a dynamical variable in our 
current model. However, by repeating the simulation for a sufficiently large range of / we are able to produce the 
curves for fixed external force / = 0.4,1.0,2.0,3.0 (from bottom to top). We assume that because there are pinning 
sites distributed along the filament, the cylinders can move through a tunnelling or hopping mechanism between sites 
and so I and hence {d±){l) will vary to minimize the free energy Thus we can infer from $ the effective interaction 
between two cylinders in a sytem where I is a dynamical variable. We choose to plot $ versus the value of {d±){l), 
the projected horizontal separation between the two cylinders, since this is the more relevant observable. One can see 
the looped phase and the extended phase, corresponding to negative and positive {d±), respectively. In the extended 
phase $ increases with the magnitude of the external force for both kinds of wrapping symmetry, as one would expect. 
The effective interaction between the cylinders does, however, depend crucially on the symmetry of wrapping. In the 
asymmetric case (Fig. [9j left pane) the effective interaction is attractive in the extended phase. In the symmetric case, 
however, (Fig. [o] right pane), we see that the effective interaction is repulsive in the extended phase but then changes 
sign in the looped phase. In the extended phase, therefore, the cylinders with symmetric (antisymmetric) wrapping 
will move towards smaller (larger) {d±); there is an effective attractive (repulsive) force between the cylinders. The 
higher the external force the bigger this effective force is between the wrapped cylinders. 



In Fig. 10 we plot the free energy E{l,{d±),ai,a2) defined in Eq. (51 1; S is a function of {d±) at fixed I 
(i.e., / = AttR in the figure). In the lefthand plot we show S for both symmetric and antisymmetric wrappings for 
/ — 0.4, /i = 50, a = 13.0 and a = 5tt/8, 3tt/8. The significance of the values chosen for a is that they are, respectively, 
^ 7r/2. For symmetric wrapping the equilibrium (minimum of S) is in the looped phase for a < 7r/2 but it is in 
the extended phase for a > 7r/2. In contrast for antisymmetric wrapping, equilibrium is in the extended phase for 
both wrapping angles although at a smaller value of {d±){l) for a < ■n/2 than for a > tt/2. In the righthand plot we 



show the projected separation {dj_) versus its conjugate variable A. From Eq. (51 1 we note that A is determined as a 
fimction of {dj_) via 

In an extended phase, {d±) > 0, we have that A > (or A < 0) corresponds to an external force applied between the 
cylinders with magnitude |A| which pushes them together (or pulls them apart). In a looped phase, (dj_) < 0, we have 
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FIG. 10: Left: the free energy as a function of the horizontal displacement projection {d±) (normalized to dmax ~ I + 2R). 
Right: the required externally applied force, A, as a function of the horizontal displacement projection {d±}. Here we consider 
the case of constant wrapping with a = 5tt/8 > 7r/2 and a = Stt/S < 7r/2 for / = 0.4, fi = 50, a = 13.0 and / = AnR. 



that A > (or A < 0) corresponds to an external force applied between the cylinders with magnitude |A| which pulls 
them apart (or pushes them together). We conclude that X{d±) > (or X{d±) < 0) corresponds to an attractive (or 
repulsive) force between the cylinders. 

These four choices for the variables {d±) and A are shown separated by the dotted lines in the plot. The equilibria 



of the lefthand plot in Fig. 10 correspond to A = in the righthand plot. Clearly, for symmetric wrapping and a 



wrapping angle a ^ 7r/2 the equilibrium value of the projected length is {dj_) ^ giving a high probability for the 
cylinders to interact. 

One should not forget here that in reality there may be other interactions between the wrapped cylinders (as in 
nucleosomes) that are not taken into account in this model: there would be direct electrostatic repulsions acting in 
real 3D space, as well as short range steric interactions when the cylinders are wrapped symmetrically, but not when 
they are wrapped antisymmetrically. It is the sum of all these complicated interactions that would need to be taken 
into account in a complete theory of cylinder wrapping. 



B. Unconstrained wrapping angles 



Next we consider the case of two cylinders with no constraints on the wrapping angles, corresponding to the second 
case (case II) as explained in the beginning of Section IV and IV B This case is more complicated from the previous 
one as it entails an additional integration with respect to the two wrapping angles. The partition function can then 
be written in the same way as before in Eq. (41 ) 



daida2 



min(//i?,Qi) 



da'i 



inin{/ / R~a'^ ,ck2) 



da'2 Z\{ai,a2,l = I' — Ra'i — Ra'2), 



(53) 



and the variables are defined in Section IV B following Eq. (41 1. We consider the free energy, now for unconstrained 
wrapping angles, given by 



(3<i>{l',X) = -\nZx{l'). 



(54) 



and its Legendre transformation 



S(/',(dj_)) =mm($(r,A)-A(di)), 



(55) 



which again correspond to a system where either the polymer length, or the average horizontal displacement, {d±), 
are fixed respectively. 
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FIG. 11: The free energy for the case of unconstrained wrapping angles. Left: the free energy is shown as a function of external 
force / for a single cylinder and for two cylinders with symmetric wrapping. The upper dashed curve is for I' = nR and the 
lower dashed curve corresponds to /' — 20nR, being effectively infinite. As one expects, for large separation and low force, the 
free energy of two cylinders is double that of a single cylinder (red curve). Right: the free energy as a function of projected 
horizontal distance d±, for / = 0.01, 1,2,3 (from bottom to top). The free energy decreases with separation (entropic effect) 
and increases with force. In both cases n = 1 and a — 1.25. 



In the lefthand pane of Fig. 11 we plot the free energy, $, of a quenched system (see remarks for case II in the 



beginning of Section IV I as a function of external force / for a single cylinder and for two cylinders with symmetric 
wrapping and various arc-length separations. For antisymmetric wrapping the results are almost identical and we do 
not show them separately. We see that as / — > and the arc-length separation /' becomes very large {{d±) — )• oo), 
the free energy of two cylinders is twice the free energy of a single cylinder, which is a good consistency check. Note 
that for large / the elastic filament can unwrap which allows the energetically favored small wrapping angle: The 
cylinders are pinned with no significant length of filament wrapped on them. This is in strong contrast to the case 
of constrained wrapping angle shown in Fig. |8] which shows a desorption transition due to the constraint imposed 
on the wrapping angle. In the righthand pane, we show the dependence of the unconstrained free energy, 4>, on the 
projected separation {dj_). Although I' is a quenched variable, we are again able to produce these curve by repeating 
the numerical calculation for a sufficiently large range of /'. The cylinders can move through a tunnelling or hopping 
mechanism between pinning sites and so V and hence {d±)(l') will vary to minimize the free energy $. As the projected 
separation between the cylinders becomes smaller, i.e. they get closer together, we get an increase in the free energy 
corresponding to effective repulsive interactions which are not due to any hard-core repulsion between the cylinders, 
but are entropically generated. This entropy stems primarily from an "entropic wrapping" effect: By limiting the 
space between the cylinders they cannot wrap in as many ways as for large separations. Such entropic wrapping 
effects should be distinguished from the usual entropic configuration effects in semi-flexible polymers. In the latter 
the number of configurations of the polymer chain changes as we restrict the position of its ends and this leads to 
entropic polymer elasticity; for entropic wrapping effect the physical picture is altogether different. 

In Fig. 12 we plot the conjugate variable A against the horizontal displacement {d±) for a system with I' = 2ttR. 
Again A can be interpreted as an external force needed to maintain a mean projected separation {d±), as described 
in Section |V A| The response of the system to the external force is monotonic and attests to the fact that a large 
repulsive force, A > 4, greater than the external applied force, / = 3 in this case, would be required to sustain a 
looped phase. The choice I' — 2ttR is relatively short and we do not include the exclusion of one cylinder by the other 
which would restrict the possible configurations especially in a looped phase. However, we treat the result here as 
an idealized case which is indicative of the possible likely configurations; in practice, the cylinders are not of infinite 
extent and so can be assumed to pass by each other more readily than the cylinders of this model. The absence of 
a looped phase in this idealized case is strong evidence that it is not likely to occur except for large A in a more 
realistic model of unconstrained wrapping. It should also be noted that the amount of computation time to explore 
all the parameter space in this case is considerable and so we are confined to investigating whether or not there is any 
significant non-trivial configuration that is likely to be realized. 

We do not plot an analogue of the lefthand pane in Fig. [TO] as it does not provide any extra information, however 
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FIG. 12: The required externally applied force, A, as a function of the horizontal displacement projection (dx) (normalized to 
dmax = I + 2i?) is shown for the case of unconstrained wrapping with / = 3, /i = 1, ct = 1.25 and I' = 2tiR. 



the free energy, S, can be easily calculated numerically in our formalism. 

From the numerical solution presented above it thus follows that the cylinder wrapping and the associated entropy 
presents yet another, apparently more important, source of polymer-mediated interactions between wrapped cylinders. 
Its source is the wrapping degrees of freedom that are constrained as the cylinders move closer together. To our 
knowledge, this source of effective interactions along a polymer chain has not before been clearly discussed in the 
literature. 



VI. SUMMARY AND CONCLUSIONS 



We have analyzed wrapping equilibria of one and two cylinders on a semi-flexible filament driven solely by the 
elastic energy of the filament and the (adhesive) energy of wrapping around the cylinders. We derived the statistical 
properties and the free energy of wrapping in the case of one and two cylinders as well as the effective interaction free 
energy between two wrapped cylinders along the elastic filament. Our calculation is based on the functional integral 
representation of the partition function for the filament and is exact, within the confines of the worm-like chain model, 
the assumed form of the wrapping potential and the limit of a 2D Eulerian plane. 

In Section [ll] we presented the generalized theory for elastic filament wrapping on one or more cylindrical cylinders 
in one dimension. In Section III wc calculated the wrapping transition for a single cylinder and showed that it is 
necessary to solve the full Schrdinger equation in order to obtain a good numerical value for the critical unwrapping 
external force. In Section |IV| we analyzed the exact solution for two cylinders pinned a fixed length apart on the 
elastic filament for both the looped phase, where the mean projected separation cZ_l is negative, and for the extended 
phase, where d± is strictly positive and there are consequently no loops. The two cases considered are of constrained 
and unconstrained wrapping, respectively. In Section [IV A| the wrapping angles ai and a2 on the respective cylinders 
are fixed: the case of constrained symmetric wrapping. This gives the most interesting results concerning the presence 
of a looped phase as discussed already by Rudnick and Bruinsma [S]- We considered two values of wrapping angles: 
oil = a2 — and ai — a2 — Stt/S. For the larger value there is a clear looped phase shown in Fig. [s] characterized by 
the average horizontal separation {d±) < 0. The loop initially increases in size as the external tension / is increased 
and eventually for sufficiently large / the system switches over to the extended phase. In contrast, for smaller value 
of ai shown in Fig. [6] there is no looped phase. In both cases the largest possible extension is close to its maximum 
possible value, dmax as determined by the arc-length of filament between the two cylinders and their radii, as one 
would expect. We conclude that a looped phase is possible for constrained symmetric wrapping and sufficiently large 
wrapping angles but is absent if the wrapping angle is too small. From Fig. [5] we see that where it does occur, the 
maximum loop size increases as the rigidity fi increases and is sustained for a range of tensions /; for fi — 10 the most 
negative val ues of {d±) are approximately for 0.2 < / < 1. 

In Section IV B| the two cylinders are pinned a distance /' apart along the contour of the elastic filament and the 
wrapping angles ai are now dynamical (annealed) variables. In contrast to the constrained case the wrapping energy 
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encoded in the dimensionless variable a plays a direct role in the values of the observables. The wrapping angle on 
a given cylinder is divided into two parts, which are the wrapping angles of the elastic filament wrapped to the left 
and to the right of the pinning point. In the two-cylinder case there are then four dynamical angle variables over 
which to sum, and this greatly increases the computer time required to carry out the calculation. We discussed in 
detail only the symmetric wrapping configuration as the antisymmetric one gives similar results. In Fig. [7] we show 
both {d±) and (a) (normalized to their maximum values) as a function of / for a = 1, 10 corresponding to small 
and large wrapping energy. For a = 1 there is a clear unwrapping transition for the two values of rigidity /i = 1, 10 
as the external tension / increases, but for larger value of a the system remains maximally wound for all values of 
/. In particular, there is no looped phase indicated by {d±) < 0; this is characteristic of the unconstrained model. 
We chose separation /' = 27ri? which is relatively short compared with the cylinder radius R and treat the results as 
an idealized case which is indicative of the possible more realistic configurations. In practice, the cylinders are not 
cylindrical or of infinite lateral extent and so can be assumed to pass by each other more readily than the cylinders 
in the present model. 

In Section [lV C| we calculate the free energy of the symmetric and antisymmetric constrained systems of two cylinders 
pinned a distance I apart with wrapping angles ai = ±a2 = tt, respectively, and compare with the free energy of 
a single pinned cylinder with a = tt. In all cases the free energy is normalized by subtracting that of the elastic 
filament with no pinned cylinders. The results show that for large I and for both the symmetric and antisymmetric 
configurations, we observe that desorption occurs for single as well as double cylinder systems simultaneously and for 
the same external tension. For small I the situation is however different. For the antisymmetric configuration the 
double cylinder system is more strongly bound than for a single cylinder and remains bound after the single cylinder 
has been desorbed. In contrast, for the symmetric configuration, as the force increases, first the double cylinder 
becomes unstable leading to a single cylinder desorption, leaving a bound single cylinder which then desorbs as the 
external tension is increased further. The conclusion is that the symmetry of the double cylinder constrained wrapping 
has a crucial effect on the desorption transition. For unconstrained wrapping there is little structure since we do not 
associate a significant binding energy with the pinning site itself and so, as the external tension increases, the cylinders 
simply unwrap but remain pinned nevertheless. The effect in the constrained case is due to the competition between 
the wrapping and the entropic contributions to the free energy as a function of external traction and separation 
between the cylinders. 

In Section |V] we investigated the induced force between two pinned cylinders by introducing a force A conjugate to 
the projected distance d±. We conclude that in the case of constrained wrapping the effective force, given by the slope 
of the curves, depends on the symmetry of configuration, being repulsive in the symmetric case and attractive in the 
antisymmetric case. The dependence of the effective force on the arc-length separation I between the cylinders follows 
closely the dependence on the projected separation between the cylinders in the direction of the external tension. In 
contrast, for the case of unconstrained wrapping we observe repulsion both for the symmetric as well as antisymmetric 
configurations. We interpret this repulsion as due to wrapping entropy that depends on the separation between the 
cylinders. This entropy differs from the usual polymer conformational entropy and one should distinguish between 
the two. The identification of the wrapping entropy presents a new concept in the analysis of the entropic effects in 
the context of polymer-particle complexes. 

We finally calculate the force A required to sustain a given mean projected separation {dj_) and derived the free 
energy {d±)). For constrained wrapping we considered both symmetric and antisymmetric configurations with 
«! = ±a2 = 37r/8, 57r/8, respectively. In Fig. [lO] we show S and A versus {d±) / dmax- It should be noted that 
\{d±_) > 0(< 0) corresponds to an intrinsic attractive (repulsive) force between the cylinders caused e.g. by charges 
on the cylinders. The behavior of the A versus {d±)/dmax curves is consistent with this interpretation. We conclude 
that for the given choice of parameters the looped phase only occurs for symmetric wrapping and ai — a2 > tt/2. 
The range of A chosen includes values where its magnitude exceeds the value of the applied external tension / = 0.4. 
The response of {d±) to A is as expected and we see a looped phase for a sufficiently repulsive intrinsic force. Other 
parameter choices can be investigated but we do not present the results here. In comparison, the case of unconstrained 
wrapping is basically featureless and the results are shown in Figs. |ll| and [T2l where only symmetric configuration is 
considered, the results for antisymmetric configuration being very similar. We note that the force between cylinders 
is repulsive in both these cases. We interpret this repulsion again as due to the wrapping entropy that depends on 
the separation between the cylinders. Since the wrapping entropy might also play an important role in the case of 
multiple wrapped cylinders and could promote very strong non-pairwise additive effects, we plan to study its effects 
very carefully in the future. Also we intend to introduce a chemical potential for exchange of the wrapped cylinders 
with a bulk phase in order to generalize the calculation of the distribution of nucleosomal core particles within the 
genomes [T^IT^ . 

A major conclusion of our work is that for constrained wrapping, where the amount of elastic filament wrapped 
around the cylinder subtends a fixed angle at the center, there are two kinds of transition that can occur as a function 
of the dimensionless external tension /, rigidity ^ and wrapping energy a. For two or more wrapped cylinders there is 
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a transition from a looped to an extended phase which is additionally affected by the direct inter-cylinder forces, and 
there are desorption transitions which are sensitive to the symmetry of the wrapping (determined by the relationship 
of the signs of wrapping angles) and also to the inter-cylinder separation. 

A second major conclusion is that for unconstrained wrapping neither the looped phase nor the desorption transition 
are likely to exist. Instead, there is an unwrapping transition where amount of filament wrapped on each cylinders 
rapidly decreases as the external tension / passes through a critical value. Correspondingly, the inter-cylinder distances 
rapidly increase from small to near maximum values within a very small tension interval. 
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Appendix A: Constrained wrapping expressed in terms of the Mathieu functions 



The path integral K obeys the Schrodinger equation 

9X(V,V',0 



-HK, (Al) 



dl 

with boundary condition 

if(V;,V^',0) = (5(7^-^')- (A2) 

This clearly means that K{il> -I- 2mT,ip',l) ^ K{tp,tp',l) as it is violated aX I = in the initial conditions. However, 
the propagator Km derived using the Mathieu functions, 5'„, 

KMi^Jj,i^',L) = ^exp(-i;„/)*„(^)«'„(7^') (A3) 

n 

has initial conditions 

ifAf(V',V'',0) =^<S(V-V''-2n7r), (A4) 

n 

and is clearly periodic. We can thus write 

XAf(V',V'',0 =I]^(^,V'' + 2n^,0- (A5) 

n 

If we take a single cylinder with fixed a (this is crucial in the argument that follows) we have 



Z{a)= j d'iJjodipidtlj3 K{il;o,tljiJi)S{a,'4ji)K{iji + a,tlj3,L ~ li - R\a\), (A6) 

where S* is a general boundary terms which is periodic in ipi. The initial integral over ipQ can clearly be taken over 
the interval [0, 2tt]. If we restrict the integrals over ipi and 1^2 to [0, 27r] and add on their integer changes by hand we 
get 

Z{a) = / dipodtpidips^'^^ K{ipo,ipi + 2mT,li)S{a,ipi + 2mT) X 

•^0 m,n 

X K{il;i + 2mi + a, V^3 + 2m7r, L - /i - R\a\), (A7) 
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where we explicitly show that we restrict the integrals to [0, 27r]. Using the fact that S is periodic for fixed a we 
finally obtain 

Z{a)^ I d^odilJidi;3'^K{il;o,tpi+2nnJi)S{a2,-ipi)K{-tpi + 2nTT + a,il;3 + 2mTT,L-h- R\a\), (A8) 

^ m,n 

and the obvious relation that 

K{ijj + 2m:, ^p' + 2rrm) = K{ij, ^p' + 2(m - n)7r). (A9) 
We then change the summation variable over to to to — n to obtain 

Z{a) = J d-tpod-ipid-ips I^^Ki-iljo,^^! + 2nTr,li)^ S{a,^i) X 

X ^^K{iPi+a,^3 + 2m7r,L-h-R\a\)^ , (AlO) 



which then gives 



^27r 

Z{a) ^ I d2pod^Jid'il}3KM{^JQ,iJi)S{a2,ipi)KM{ipi + c(,ip3,L - h - R\a\), (All) 



which is the desired result expressed in terms of periodic Mathieu functions. The proof above for fixed wrapping 
angles can easily be extended to several cylinders with fixed wrapping angles. 

Appendix B: Horizontal distance between the cylinders 



From Eq. ( 40 1 we can write for the horizontal distance between the cylinders 

d± = (rfi) + i?sgn(ai)(sin(V'2)) - i?sgn(a2)(sin(V'3)), (Bl) 

where 

(rfi) = 777 — ^ n I I dil;idip2'^a{tpi)CaJipi)'S^exp{~e„,s/R)exp{-en{l~s)/R)x 

Z{ai,a2j)J J ^ 

X y"d(sinV's)^'™(7/'^)*™(V'i +ai)^'„(V's)*n(V'2)C„3(V'2)^'o(^/'2 + a2), (B2) 

where a state has been inserted, with angle ips at the point s of the first cylinder. We propagate the solution up to 
this point, calculate the horizontal projection, and then propagate to the remaining cylinder. We now find that 

(rfi) = 777 — A / dipidtp2'ifo{ilJi) ■ Ca,{->Pi)Y\ / ds exp(-e„s/i?) X 

Z{ai,a2,l)J 

X exp(-e„(; - s)/R)Dmn{i^i + "l, V'2)Ca2(V'2)*o(V'2 + ^2), (B3) 

where 

-Dmn(V'l + «!, ^"2) = / C^(sinV's)P,„(-0i +ai,-0s)P„('0s,V'2), (B4) 



and Ca is defined in Eq. (37). By noting that after the angular integration only the exponential dependence on s 



remains, we integrate s over the range < s < / to finally get 

^Ze-^"'/^7D„„(^i+a2,^2)- 



(C?i> = 77^^ ^ / rfV'ld^2*o(V'l)ai(V'l) 



Z{ai,a2,l) 



■ 2^ Aim (-01 + "2,-02) 



Ca,(V^2)*0(V'2+a2). (B5) 
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Note that because L ^ I, the effect of the filament lengths outside the cylinder region is encoded in the ground state 
factors \['o(V'i) and ^'o(i/'2 + 0:2)- The average distance {d±) is easily computed since the various components can be 
pre-computed. In a typical computation we pre-compute the lowest 20-100 Mathieu eigenfunctions and eigenvalues 
by recasting the Schrodinger equation in Eq. ( 29 ) as a matrix eigenvalue problem by discretizing the angle coordinate 
on the range [0,27r], and using proprietory NAG routines j44| . The computation of the partition function is then 
straightforward and can be done with modest computing resources. 
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